Universal Magnetic Fluctuations with a Field Induced Length Scale 
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We calculate the probability density function for the order parameter fluctuations in the low tem- 
perature phase of the 2D-XY model of magnetism near the line of critical points. A finite cor- 
relation length, ^, is introduced with a small magnetic field, h, and an expression for ^(/i) is de- 
veloped by treating non-linear contributions to the field energy using a Hartree approximation. 
We find analytically a series of universal non-Gaussian distributions of the finite size scaling form 
P(m,L,0 - L'^^^PLimL''^'" ,i/L) and present a function of the form P{x) ~ (exp(2: - exp(x)))"(''' 
that gives the PDF to an excellent approximation. We propose a{h) as an indirect measure of the 
length scale of correlations in a wide range of complex systems. 
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I. INTRODUCTION 

There has recently been considerable interest in the 
fluctuations of a spatially averaged quantity in systems 
with correlations over a macroscopic length scale ^ . 
The most accessible example from both an experimen- 
tal and a theoretical point of view is a critical system, 
where the divergence of the correlation length ^ is in- 
terrupted by the system size L. Here, the breakdown 
of Landau theory and the occurrence of a non-analytic 
fixed point for the free energy in a renormalization group 
flow, tells one that the fluctuations will not be Gaussian. 
Even so, extending the scaling hypothesis to include the 
macroscopic scale L, one can deduce that the probability 
density function (PDF) for order parameter fluctuations 
should have universal properties. Further, the PDF must 
be a homogeneous function of to, L and ^ of the following 
form 

Pim.L.O-L^'" PdmLP'-.^/L). (1) 

Here, adopting the language of a magnetic phase tran- 
sition, V and (3 are the usual critical exponents relating 
to the divergence of the correlation length and the sin- 
gularity in the magnetization, m. The scaling hypothesis 
therefore predicts fluctuations of a universal form, inde- 
pendently of system size, for constant ratio ^/L. 

We study the low temperature phase of the 2D-XY 
model, defined by the Hamiltonian 

H = -J^cos{e,-ej) -h^cosOi. (2) 

The exchange interaction and magnetic field are of 
strength J and h and the angle 9i gives the orientation of 
a classical spin vector of unit length, confined to a plane. 
We define the magnetization for a single configuration 

m=l ^ cos(0,-0), (3) 



where 9 = tan"^ sinfl^/ cos 6*^) is the instanta- 
neous magnetization direction. 

This is perhaps the simplest nontrivial system in which 
one can study critical phenomena. At low-temperature 
and in zero field there is a line of critical points, sep- 
arated from the high temperature paramagnetic phase 
by the Kosterlitz-Thouless-Berezinskii phase transition. 
The physics of this low temperature phase is perfectly 
captured by a harmonic, or spin wave Hamiltonian. That 
is, one can, without loss of generality, develop the cosine 
interaction to order {9i — 9jY and neglect the periodicity 
of 9i. This Hamiltonian is diagonal in reciprocal space 
and can be solved straightforwardly. As a result, all crit- 
ical phenomena can be calculated microscopically from 
Gaussian integration, without the need for either the 
scaling hypothesis, or the renormalization group. Along 
the line of critical points the exponents [3 and v are not 
individually defined, but there ratio is and the system 
has a single independent exponent rj = 2f3/v = T/2'kJ. 

We have previously been interested in the zero field, 
or strongly correlated regime where the divergence of ^ is 
completely removed by the system size L and Pl becomes 
a function of a single variable mL^/'^ . We have found 
that Pl, when plotted as a function of /x = (m— {m))/a, 
is a universal function, not only of system size, but also of 
temperature and therefore of critical exponent 77. Here, 
(to) is the mean and a the standard deviation of the 
distribution. This rather surprising result gives weight 
to our conjecture ||] that the critical fluctuations of sys- 
tems in certain universality classes are captured, at least 
qualitatively by the fluctuations of the 2D-XY model. 

In this paper we generalize our previous results by 
introducing a second length into the problem with the 
aid of a magnetic field. The field breaks the symme- 
try moving the system into an ordered magnetic state 
with finite correlation length ^. However, taking a van 
Hove type thermodynamic limit, with the ratio ^/ L con- 
stant, should lead to a family of limit functions, all with 
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divergent correlation length, varying in form from the 
anisotropic limit (see figure 1) to a Gaussian function, 
as the ratio ^/L falls to zero. In the next section we 
develop a starting Hamiltonian that satisfies the require- 
ments of the scaling hypothesis, using a self-consistent 
Hartree approximation. In section 3 we give theoretical 
results for the PDF for finite field and compare our re- 
sults with Monte Carlo simulation. In section 4 we fit the 
curve with a generalized form of Gumbel's first asymp- 
tote from extremal statistics. 



II. HARTREE APPROXIMATION AND THE 
HYPERSCALING RELATION 

Expanding the field energy in small angles, in the same 
way as the exchange term and Fourier transforming gives 
a convenient starting Hamiltonian 



-E 



h 

lq + -j 



(4) 



where (pq — Re l/vW 0^ exp(— iff.fi) is the real 
part of the Fourier transformed spin variable. 7^ — 
4—2 cos{qx) — 1 cos{qy) and the sum J2q^o '^^'-^^ ^^'^ Bril- 
louin zone for a square lattice with periodic boundaries, 
with q taking on discrete values qx — {2Tr/ L)nx, qy = 
{2tt /L)ny, nx,ny = 0,2...y/lN). Here and throughout 
the paper we have set the nearest neighbour distance on 
the lattice equal to unity. Expanding 7^ for small q we 
can write the Green's function propagator 



G{q) 



(5) 



the magnetic field indeed introduces a length scale 

However, this naive starting point needs some devel- 
opment before proceeding with the calculation as, as it 
stands it does not satisfy the well known hyperscaling 
relation. To see this, consider the following scaling argu- 
ment I p^ : at the critical temperature but in finite field, 
the thermally averaged magnetization can be expressed 
in terms of both ^ and h 



-13/1^ 



J 



1/5 



(6) 



where S is the usual critical exponent for the response 
in finite field. Putting the expression for ^ in (|^) 
leads to a relation between the exponents: 6 — 2iy/P, in 
disagreement with the hyperscaling relation S+1 = dv/ f3, 
which should be vahd for the 2D-XY model Hi. The 
error comes from the development of the field term in 
small angles. Even at low temperature, when the nearest 
neighbour differences 9i — 9j are small, the deviations 
of 6i from the fixed field direction are divergent in the 
thermodynamic limit. The development of the field term 



in small angles is therefore invalid. This problem can be 
dealt with in the low temperature phase, in the absence of 
vortices, using the Hartree approximation introduced by 
Pokrovsky and Uimin [p^ . Expanding cos{9i) in powers 
of 6*1 we make a mean field decoupling 



2p 



r 



.2p-2\ 



(7) 



where C„ 



(ML 



is a binomial counting factor. As the 



~ 2(2p-2)! 

underlying Hamiltonian (^ is quadratic we can reduce 
{9^^~'^) using Wick's theorem and after some resumma- 
tion we eventually find 



cos 9i Ki 1 — (m) 



(8) 



The field term in the Hamiltonian (^ then becomes 
hJ2cosi9,)^Nh-lKff{T)J24'l 



iT-effiT) = {m)h. 



(9) 



Using the scaling relation (m) ^ (h/ J)'^ the effective field 

heff — > and the above scaling argument correctly 

yields the hyperscaling scaling relationship defined above. 
Note however, that the scaling argument is valid in the 
thermodynamic limit where ^/L <IC 1 and the influence of 
the finite system size is negligible. In the crossover region 
that interests us, with £,/L 0(1) one cannot make this 
substitution and in general one must explicitly work with 
expression (^. 

The point of principle that poses the problem for the 
hyperscaling relation is that making the substitution 
cos{9i) w 1 — {l/2)9f results in an order parameter conju- 
gate to the field, (m) = (1- (1/2)6*2) ^ l -T/Stt Jlog(iV) 
which diverges with system size, for any finite tempera- 
ture. In order for the hyperscaling relation to hold, (m) 
must be a correctly defined intensive variable. For this 
to be so, the higher order terms in the expansion of cos^i 
must be retained, at least within the level of approxima- 
tion shown here. The calculation of Berezinskii p^ ] is 
consistent with this thermodynamic argument. 

Still in the absence of vortices, an effective coupling 
constant Je// can be calculated in a similar manner. In 
zero field we did not need to take account of this, as 
Pl for a quadratic Hamiltonian is independent of tem- 
perature throughout the low temperature regime. How- 
ever, as the correlation length depends on both the field 
and the coupling constant, we now have to calculate it 
if we are to have good agreement with numerical data. 
Expanding cos{9i — 9j) in powers of the discrete differ- 
ence operator V9i [|l^ and again using the decoupling 
(0) for (V6>,:)2p, we arrive at a self-consistent expression 
for J,ff{T) (l|l| 



Je// = Jexp(-^) 



(10) 
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The Green's function we finally use for the calculation 
of the PDF is therefore of the form (||) but with the 
correlation length given in terms of the self-consistent 
effective field and coupling constant 



'J. 



eff 



(11) 



The variation of ^ with temperature, for fixed field is 
quite small throughout the range of fields that interests 
us. In Figure 1 we show {£,{0) / £,{T))'^ as a function of 
temperature for three different field strengths. Even for 
h/ J = 0.5, there is only a 10% variation, up to a tem- 
perature T/J = 0.7, above which the Hartree approxi- 
mation breaks down. As shown below, the temperature 
dependence of the resulting distribution function is even 
weaker than that for £/L and for practical purposes it 
can be considered as temperature independent. 
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FIG. 1. {£.{0)/£.{T))'^ , as given by the Hartree approxima- 
tions for, from bottom to top, h/ J = 0.01, 0.05 and 0.5 (equa- 
tions (|) and (|l^) as a function of T/J. 

The parameters hef / and ^ /L can be found in table 1 
for a system of size L = 32 at temperature T/J = 0.7 
and for field strength between h = 0.001 and h = 0.5. 
For finite field the ratio ^/L varies from ^/L = 0.99 to 
C/L = 0.043. 
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-0.89 
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0.001 


0.987 


1.611 


-0.88 


0.005 


0.005 


0.4407 


1.7416 


-0.844 


0.01 


0.0101 


0.3111 


1.903 


-0.801 


0.05 


0.0512 


0.1381 


3.359 


-0.583 


0.1 


0.1034 


0.098 


5.333 


-0.463 


0.5 


0.523 


0.0429 


21.82 


-0.23 



TABLE 1 Variation of /le//, £,/L, a and 7 with field h for 
L = 32 



THE PROBABILITY DENSITY FUNCTION 
IN FINITE FIELD 



III. 



We have previously developed ||11|,|1^ the following ex- 
pression for the PDF 



(12) 



q#0 



q#0 



1 -i 



2 G(q) 



92 N 



where gk = l/N^ J2qGiq) and fj, = {m - {m))/a. The 
PDF with a finite correlation length is calculated from 
the same expression by inserting the the modified Green's 
function (||) . The equivalence of equations (|l|) and ( p^ ) 
is a result of the hyperscaling relation ||ll| and the func- 
tional dependence predicted by the scaling hypothesis 
(|l|) comes directly from dimensional analysis of equa- 
tion (|l2|). We note that the calculation can easily be 
extended to explore the non-Gaussian but non-critical 
behaviour in all dimension less than four ||ll[]. Summing 
over the Brillouin zone for a large but finite system and 
performing a numerical Fourier transform we generate 
the data shown in figures 2 and 3. Data is shown in fig- 
ure 2, for h/J = 0.05 {£,/L = 0.138) and in Figure 3 for 
h/J = 0.01, 0.05 and 0.5 corresponding to £/L = 0.311, 
0.138 and 0.043. It is compared with results from Monte 
Carlo simulation for a system of size L = 32. In each case, 
theoretical and numerical data are shown for three tem- 
peratures T/J = 0.1, 0.4, 0.7. As in the zero field case, 
which can be considered as the extreme non-Gaussian 
limit for such a system, the PDFs are characterized by 
an exponential tail for fluctuations below the mean and 
a double exponential above the mean ||Tl|| . Applying the 
field reduces the asymmetry and in large field the data 
approach a Gaussian distribution. 

Agreement between the theoretical calculation and the 
Monte Carlo simulation is generally extremely good, in- 
dicating that the Hartree approximations are accurate. 
In Figure 2 all sets of data collapse, within numerical 
error, onto a single curve independently of temperature. 
This is the case for all field values chosen. When plotted 
on a logarithmic scale, temperature dependence is still 
not observable, but a difference between theoretical and 
numerical values can be observed along the exponential 
tail, for probability densities smaller than 10^"*. The dis- 
crepancy appears largest for fields around h/J ^ 0.05. 
This must indicate the limit of the ability of the Hartree 
approximation in dealing with the fluctuations. For very 
small field, its effect is small and so errors are negligible, 
while for larger fields critical fluctuations are smaller and 
one can imagine that the Hartree approximation becomes 
quantitatively very accurate. Only in the intermediate 
field range of h/J ^ 0.05, the combination of these two 
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effects is sufficient to give the small deviation from the 
Monte Carlo data. 

0.5 I > > > > 1 




1^ 

FIG. 2. Monte Carlo data for Pl{i-i.) for magnetic field 
h/J = 0.05 and for T/J = 0.1,0.4 and 0.7, using the Hamil- 
tonian (^. The lines are for data generated from equation 
@ (dotted) and for the function ^ (dashed). 

In all cases, both numerical and theoretical, the inde- 
pendence of the results on temperature is quite remark- 
able and it leads one to suggest that, as in the case of zero 
field, the distribution is truly temperature independent 
throughout the range of temperature and system size for 
which the excitation of vortex pairs can be neglected. 
This point requires further study, but it is already clear 
that, from a pragmatic point of view, temperature de- 
pendence is not an observable phenomena. 

Note that with the definition (^) we explicitly study the 
longitudinal magnetization, irrespective of its direction 
in space. This is the quantity that becomes critical at a 
phase transition in a system with continuous symmetry. 
The introduction of the magnetic field breaks the orien- 
tational symmetry in the thermodynamic limit and the 
variable conjugate to it is the projection of m along the 
field direction. The fluctuations in these two quantities 
are different for small field and in the crossover region, as 
h ^ 0, the latter quantity becomes ill-defined The 
two quantities become indistinguishable for ^/L ~ O.I. 

The skewness 7 = (/x"^), which parameterizes the asym- 
metry of the curve, varies from —0.89 in the extreme 
non-Gaussian limit to —0.23 for a ratio £,/L = 0.043 and 
indeed goes smoothly to zero for £^/L 0. Numerical 
values can also be found in table I. 
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FIG. 3. Monte Carlo data for Pl{i-i) for magnetic field (a) 
h/J = 0.01, (b) h/J = 0.05, (c) h/J = 0.5 for T/J = 0.1, 0.4 
and 0.7. The lines are for data generated from equation ( |l^ ) 
(full) and for the function (|l^) (dashed). Theoretical curves 
for different temperatures are superposed. 

The origin of the skewness becomes clear if one con- 
siders the contribution made by the normal modes. The 
magnetization can be written 
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(13) 



9)0 



To leading order in (j)q m therefore consists of a sum over 
variables mq = {4>^/2N) which, within the spin wave 
approximation, are statistically independent, with gen- 
erating function 



P(mq) 



^:^m-i/2e-/3^=//^(?' . (14) 
Ait ^ 



In zero field the mean amplitudes (mq) vary from a mi- 
croscopic value, 0{1/N) for modes on the zone edge 
through to a value of 0(1) for the long wavelength modes 
at the zone center and the dispersion in the contribu- 
tions is divergent in the thermodynamic limit. In two- 
dimensions the density of states is linear in q which is 
just what is required to engage the entire zone |ll| . 

Violation of the central limit theorem therefore arises 
because the individual elements, although statistically in- 
dependent are not individually negligible. The modes of 
divergent amplitude near the Brillouin zone center are 
responsible for the anisotropy although all parts of the 
zone are required for a detailed reconstruction of P{m). 
Introduction of the length scale ^ removes the divergence 
for g and re-establishes the criterion that the statis- 
tically independent elements are individually negligible. 
In the limit that ^/L ^ the distribution becomes Gaus- 
sian. If the thermodynamic limit is taken while keeping 
the ratio ^/L constant the amplitudes remains divergent, 
but the contribution from the zone center becomes pro- 
gressively less, as the ratio ^/L is reduced and the skew- 
ness falls to zero. 



IV. FITTING WITH A GENERALIZED GUMBEL 
FUNCTION 



In refs ||ll|,|l^ we have compared the functional form 
of Pl{i^) with a series of standard expressions. Although 
none are exact solutions, they all give good fits to the 
data and therefore offer very useful analytical functions 
as well as giving some insight into the physical processes 
responsible for the asymmetric PDF. Here, we only pur- 
sue one of these, the generalized Gumbel function for the 
statistics of extremes [ITsi 



(15) 



which gives the a*'* largest or smallest values of a set of 
N random numbers, z^, in the limit that N oo. For 
the smallest values, /iz = (z — {z))/az- The constants 
w, h and s depend on a through the three conditions of 
normalization, (/i^) ~ and (/i^) = 1 and one finds 



a"'aa 



T{a) 



I 1 d^T{a) 
T{a) da? 

1 

log(a) 



1 ar(a) 



r(«) 

dT{ay 



da 



r(a) da 



(16) 



The function therefore has only one parameter, which is 
calculated by comparing the Fourier transform of ( [l5| ) 
with <I'(x) of equation ( p^ ) . Using the notation of equa- 
tion dl3) we have, for the Gumbel function 



where ^{z) is the digamma function r'(z)/r(z). For $ 
we have: 



X 



(18) 



94 4 , ■ 2\/2ff5 5 , 



from which it follows that the constant a is implicitly 
given by 



3/2 



4"(a)3/2 32 



(19) 



For zero field the solution is a « 7i'/2, rather than an 
integer value, showing Pirn) is not simply an extreme 
value distribution. That this solution is an approxima- 
tion, though be it a good one, can be seen by comparing 
the ratio of higher order terms in the two expansions ( |l7| ) 
and (|l^). These diverge slowly from unity |jll|. Solv- 
ing equation (^) for finite field gives a increasing with 
h. The subsequent curves are superimposed in Figures 2 
and 3 where one can see that the fitting function repro- 
duces the results of the theoretical calculation to a good 
approximation. For small field a very small difference in 
the slopes of the exponential tails can be detected. How- 
ever, this disappears with increasing h and the fitting 
function can be regarded as an excellent working tool for 
describing the data. In Figure 4 we illustrate the evolu- 
tion of the distribution from the anisotropic limit to the 
uncorrelated Gaussian limit as a function of field using 
equation (p^. The values of a{h) are shown in table 1. In 
terms of extremal statistics, the evolution of a(h) means 
that we are describing the PDF of less and less extreme 
values, which becomes more and more normal. 
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FIG. 4. Evolution of Pl(m) from the function (|l|) with 
h/J. Moving from left to right for fluctuations below the 
mean, the data is for h/J = 0, 0.01, 0.05, 0.1, 0.5, oo 



For strong field one can solve equation(|19D analytically. 
Evaluating g2 and 53 using a continuum approximation 
and using Stirling's formula, Inr(a) « a In a— a, one finds 



(20) 



This simple expression gives a to a good approximation 
even outside the range of a values for which Sterling's 
formula is valid and reproduces our previous result for 
h = 0. It also allows one to see that a has a contribution 
coming from fluctuations within a correlated domain and 
a contribution coming from the fact that the system, with 
finite ^, can be divided into a number, N^ff = (L/^)^, of 
statistically independent domains. 



V. CONCLUSION 

In conclusion, we have made a microscopic calcula- 
tion of the generalized scaling function Pi(mi^/'', ^/L) 
for order parameter fluctuations near the line of crit- 
ical points of the low temperature phase of the two- 
dimensional XY model. A Hartree approximation is used 
to treat the non-linear corrections to a quadratic Hamil- 
tonian. The approximation is necessary to ensure that 
the hyperscaling relation between critical exponents is 
satisfied. We show that the hyperscaling relation is a 
consequence of the non-linearity necessary to ensure the 
correct system size dependence of the order parameter, 
conjugate to the applied magnetic field. This is a require- 
ment of thermodynamics, rather than a general require- 
ment for the observation of non-Gaussian fluctuations 
for global quantities. Indeed, observation of hyperscal- 
ing in non-thermodynamic systems p9[ , could be taken 
as an indication that an equivalent phenomenology ex- 
ists. For fixed magnetic field the correlation length is 



modified only slightly by thermal fluctuations and this 
manifests itself in the function Pj^ which is essentially in- 
dependent of temperature. The fact that the exponent 
f3/i/ = AnT/J is small and and that S = 8ttT/J - 1 is 
large may be important for this observation. More work 
is required to clarify this point. 

Finally, we propose that our fitting parameter a{h) 
could be used as an experimental tool to estimate the 
correlation length scale in other correlated systems. We 
have previously made an empirical observation |^,|l2| that 
the fluctuation of global measures in other correlated sys- 
tems, both in equilibrium and out of equilibrium are very 
similar to those of the magnetization of the 2D-XY model 
in zero field. We have proposed that these observations 
illustrate at least qualitatively, universal features for cor- 
related systems from different universality classes. This 
idea, as well as alternative interpretations , could be 
tested for example, in an enclosed turbulent flow using 
the experimental set up described in , by varying the 
ratio of the power injection length scale to the enclosure 
length scale. 
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